Landau Level Collapse in Gated Graphene Structures 
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We describe a new regime of magnetotransport in two dimensional electron systems in the pres- 
ence of a narrow potential barrier imposed by external gates. In such systems, the Landau level 
states, confined to the barrier region in strong magnetic fields, undergo a deconfinement transi- 
tion as the field is lowered. We present transport measurments showing Shubnikov-de Haas (SdH) 
oscillations which, in the unipolar regime, abruptly disappear when the strength of the magnetic 
field is reduced below a certain critical value. This behavior is explained by a semiclassical analysis 
of the transformation of closed cyclotron orbits into open, deconfined trajectories. Comparison to 
SdH-type resonances in the local density of states is presented. 
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Electron cyclotron motion constrained by crystal 
boundaries displays many interesting phenomena, such 
as skipping orbits and electron focusing, which have 
yielded a wealth of information on scattering mecha- 
nisms in solids [l|, Since the 1980s, many groups 
have investigated electron transport in semiconducting 
two-dimensional electron systems (2DES), where gate- 
induced spatially varying electric fields can be used to 
alter cyclotron motion. A variety of interesting phenom- 
ena were explored in these systems, including quenching 
of the quantum Hall effect [1, 0] , Weiss oscillations due 
to commensurability between cyclotron orbits and a pe- 
riodic grating Q, pinball-like dynamics in 2D arrays of 
scatterers Q, and coherent electron focusing 0.. 

The experimental realization of graphene [8|, a new 
high-mobility electron system, affords new opportunities 
to explore effects that were previously inaccessible. In 
particular, attempts to induce sharp potential barriers in 
III-V semiconductor quantum well structures have been 
limited by the depth at which the 2DES is buried — 
typically about lOOnm below the surface [§]. In contrast, 
electronic states in graphene, a truly two-dimensional 
material, are fully exposed and thus allow for potential 
modulation on ^ 10 nm length scales using small local 
gates and thin dielectric layers l3"13|- Significantly, 
such length scales can be comparable to the magnetic 
length Ib = yhcfeB, which characterizes electronic 
states in quantizing magnetic fields. In this Letter we 
focus on one such phenomenon, the transformation of 
the discrete Landau level spectrum to a continuum of 
extended states in the presence of a static electric field. 

The behavior which will be of interest for us is illus- 
trated by a toy model involving the Landau levels of a 
massive charged particle in the presence of an inverted 
parabolic potential U{x) = —ax^. Competition between 
the repulsive potential and magnetic confinement gives 
rise to a modified harmonic oscillator spectrum 
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FIG. 1: (a) Differentiated conductance, dG/dVtg, of a nar- 
row top gate graphene device, pictured in (c). Fabry- Perot 
(FP) oscillations appear in the presence of confining pn junc- 
tions, (b) dG/dVtg as a function of B and Kg. Shubnikov-de 
Haas (SdH) oscillations are observed at high B. The fan- 
like SdH pattern is altered by the barrier: in the pp'p region 
it curves, weakens, and is washed out at fields \B\ < Be, 
Eq.Q, while in the pnp region a crossover to FP oscilla- 
tions occurs. Data shown correspond to Vbg = —70V [dashed 
line in (a)], (c) Top gated graphene device micrograph and 
schematic; top gate width is ~ 16nm. (d) Local density of 
states (DOS) in the middle of parabolic barrier, obtained nu- 
merically from Egs. fllOp .d^. The derivative with respect to 
energy, dN/de, which corresponds to the measured quantity 
dG/dVtg, is shown. Dashed parabola marks the critical field, 
Eq.®. Oscillations in the DOS modulate the rate of scatter- 
ing by disorder, resulting in the SdH effect [l^ . 



for B > Be, where m is the particle mass, Py is the y 



component of momentum, and Be — V2ma/e is the crit- 
ical magnetic field strength. For strong magnetic field, 
B > Be, the spectrum consists of discrete (but disper- 
sive) energy bands indexed by an integer n, whereas for 
B < Be the spectrum is continuous even for fixed Py. 
This behavior can be understood quasiclassically in terms 
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of transformation of cfosed cyclotron orbits into open or- 
bits, which occurs when the Lorentz force is overwhelmed 
by the repulsive barrier potential. 

For massless Dirac charge carriers in single-layer 
graphene, Landau levels subject to a linear potential 
U{x) = —eEx manifest an analogous collapse of the dis- 
crete spectrum below the critical field Be = E /vp [15]: 



EniPy) = ±vp^2nheB (l 
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where n — 0,1,2... and /3 — E/vpB. The transition 
a,t B = Be can be linked to the classical dynamics of a 
massless particle, characterized by closed orbits at B > 
Be and open trajectories at B < Be [l6| . 

A simple but intuitive picture of the spectrum ([2]) can 
be obtained from the Bohr-Sommerfeld quantization con- 
dition 



Px{x)dx ~ Trh{n + 1/2 — 7), 



(3) 



Pccix) = ^ie-U{x))yvj,-{py-eBx)^, (4) 

where xi and X2 are the turning points, Py is the con- 
served momentum parallel to the barrier, and the Berry 
phase contribution 7 is 1/2 for Dirac fermions. For lin- 
ear U{x), this gives the Landau level spectrum ([2]) for 
B > Be- As B approaches Be, one of the turning points 
moves to infinity, indicating a transformation of closed 
orbits into open trajectories. 

To realize the collapse of Landau levels in an electron 
system, several conditions must be met. First, it must be 
possible to create a potential barrier that is steep on the 
scale of the cyclotron orbit radius. Second, the system 
must be ballistic on this length scale, in order to sup- 
press the broadening of Landau levels due to disorder. 
Graphene, which is a truly two-dimensional material with 
high electron mobility, fulfills both conditions. Crucially, 
as demonstrated by the recent observation of Fabry- Perot 
(FP) oscillations in gated graphene structures [l3], trans- 
port can remain ballistic even in the presence of a gate- 
induced barrier. Thus graphene is an ideal system for 
studying the Landau level collapse. 

Transport data taken from a locally gated device 
similar to that described in Ref.[l3l are shown in 
FiglH Graphene was prepared via mechanical exfoli- 
ation and contacted using electron beam lithography 
before being coated with a 7/10 nm thick hydrogen 
silsesquioxane/Hf02 dielectric layer. Narrow (^-^16 nm) 
palladium top gates were then deposited, and the electri- 
cal resistance measured at 1.6 Kelvin. Due to the simi- 
larity between the top gate width and its distance from 
the graphene 2DES, the top gate-induced density is well 
approximated by that of a thin wire above the sample 
plane and an infinite conducting plane below it: 
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where Ctg(bg) and VJ;g(bg) are the top (bottom) gate ca- 
pacitance and applied voltage and w is the distance be- 
tween the wire and the top gate. For such a thin gate, 
the modulation of the total resistance due to the gate is 
small; in order to subtract the series resistances of the 
graphene leads, the numerical derivative of the conduc- 
tance with respect to the top gate voltage, dG/dVtg, was 
analyzed. 

At zero magnetic field (Figllji), dG/dVtg shows distinct 
behavior in four regions in the (Vbg, Vtg) plane, corre- 
sponding to pp'p, pnp, npn, and nn'n doping, where n 
(p) refers to negative (positive) charge density and ' in- 
dicates different density. The appearance of FP interfer- 
ence fringes when the polarity of charge carriers in the 
locally gated region and graphene leads have opposite 
signs indicates that the mean free path is comparable to 
the barrier width, Zmf ~ w. 

In high magnetic field, in both the bipolar and unipo- 
lar regimes, we observe a fan of SdH resonances corre- 
sponding to Landau levels (see FigHJ)). At lower fields, 
different behavior is observed depending on the polarity 
under the gate. In the bipolar regime, as B is lowered, the 
SdH resonances smoothly evolve into FP resonances. The 
half-period shift, clearly visible in the data at S 1 T, 
is a hallmark of Klein scattering 
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regime, the SdH resonances bend, becoming more hori- 
zontal at lower field. The oscillations first begin to lose 
contrast, and then completely disappear below Be ~ 5 T. 

The connection between this behavior and Landau 
level collapse is exhibited most clearly by a semiclassical 
analysis. In such analysis, SdH oscillations arise from the 
interference contribution to the density of states at the 
Fermi level due to closed trajectories; as a result, Eq.Q 
with e = ep and Py = Q gives a good estimate for the po- 
sitions of the transport resonances (see also FigHJl). For 
a generic barrier potential, the quantization condition 
can be written directly in terms of experimental con- 
trol parameters. Using the Thomas-Fermi approxima- 
tion, and ignoring the effects of 'quantum capacitance' 
and nonlinear screening [l7| . we define the position- 
dependent Fermi momentum kp{x) = y^4:TTp{x)/ g, where 
g — 4: is the spin/ valley degeneracy. Substituting the re- 
lation e — U{x) = hvpkp{x) into Eq.(|3]), we obtain 



^A'nh?p{x)/g - {py 



eBx)'^dx — Trhn. 



(6) 



CbgKg 



(5) 



We see that the resonance condition is controlled solely 
by the density profile p{x), Eq.(l5|), and is insensitive 
to the specifics of the single-particle Hamiltonian up to 
slight modifications due to degree of pseudospin degen- 
eracy and Berry phase. In particular, the quantization 
condition assumes the same form, apart from the Berry 
phase contribution, for massless Dirac particles (mono- 
layer graphene) and massive particles with g — 4: (bilayer 
graphene) ; the spectrum would be only trivially modified 
for GaAs quantum wells {g = 2 and 7 = 0). 
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A simple estimate for the critical field can be ob- 
tained by comparing the curvature of p(x) at x = with 
the B'^x^ term in Eq.®. Using the device parameters 
Cbg = 115aF//^m^, Vj^g = —70 V, and w ~ 50 nm, we 

1 /2 

find i?c = {h/ew) {nChgVhg/e) ' w 5.2 T, which matches 
the data quite well (see FiglTlD). Here we used the rela- 
tion CbgVbg -I- C'tg^tg — 0, which defines the boundary 
at which polarity reversal occurs (white dashed line in 
FiglBi). 

The actual density profile p{x) is nonparabolic, flatten- 
ing out on a length scale 2?« « 100 nm. However, since 2w 
far exceeds the magnetic length for the fields of interest 
{B > IT), this flattening does not signiflcantly impact 
our discussion of the collapse phenomenon. Indeed, al- 
though the states realized at subcritical magnetic fields 
are not truly deconfined due to cyclotron motion in the 
region outside the top gated region (TGR), |a;| > w, the 
corresponding orbits are very long. For such states, the 
particle traverses the TGR, makes a partial cyclotron or- 
bit outside of the TGR, and finally crosses the TGR again 
to close the orbit (Fig. 2a). The net orbit length is a few 
w, which is much greater than the orbit size at strong 
fields (i.e. a few magnetic lengths). The contribution 
of long orbits to SdH oscillations will be suppressed due 
to spatial inhomogeneity and disorder scattering; hence 
the distinction between confined and deconfined orbits 
remains sharp despite the fiattening of the potential. 

To estimate the critical field as a function of exper- 
imental control parameters T4g and Vbg for the density 
profile ([5|), we consider an equation for the turning points. 
Setting py = 0, we have kkpix) = zLeBx. Solving this 
equation and equating the result to barrier half-width, 
Xi(2) = iw, we obtain 



B, 



= {h/ew)^{2n/eg){2C^gV^g + CtgKg). (7) 



The observed behavior of Be is well described by this 
result (red line in FiglTJj). 

To further explore the effects of curvature, below we 
analyze the inverted parabola model, U{x) — —ax^. A 
simple estimate of the collapse threshold can be obtained 
by considering balance between the Lorentz force and the 
force due to the electric field, vpB = —dU/dx. This 
condition is satisfied for a particle moving parallel to the 
barrier with x = £ = evpB / (2a). Thus we find an energy- 
dependent critical field. 



Bc{s) = {2/evFW-ae, 



(8) 



which increases with detuning from neutrality, as in ex- 
periment. 

A microscopic model of collapse is provided by the 
Hamiltonian 

( U{x) vfP- \ -fc ^ I V n ^ 

^^[vFP+ U{x) j' P±-'^^Yx^<Py-eBx\ 

(9) 



{fevlaY/\ X, 



where Py is the conserved canonical momentum compo- 
nent parallel to the barrier. We nondimensionalize the 
problem using "natural units" 

\ a J e \VFii } 

For each value of py and magnetic field B, we represent 
the Hamiltonian as an M x M matrix defined on a grid in 
position space, with periodic boundary conditions. We 
use the eigenvalues and eigenstates obtained from diag- 
onalization to evaluate the local density of states (DOS) 
in the middle of the barrier, 



M 



N{e) ^ I 'k}L^^ (IV'n,p„(a; = 0)|2) 



2^ {e-enY+l^ ' 



(10) 



with Landau level broadening incorporated through the 
Lorentzian width 7 = 0.2e*. In our simulation, a system 
of size L = 15x* discretized with M = 1500 points was 
used. Averaging with a gaussian weight was used to sup- 
press the effect of spurious states arising due to a vector 
potential jump at the boundary, 

(|V'„,p,(x = 0)n ^ / dx'e-^''/^'^'\4>n,pM)?^ (11) 



with tj w X*. Oscillations in the density of states ([TU]) 
modulate the rate of electron scattering by disorder, and 
thus show up in transport quantities measured as a func- 
tion of experimental control parameters, as in the canon- 
ical SdH effect {H. 

The resulting local DOS, shown in FiglTJi, exhibits os- 
cillations which track Landau levels at high B. At lower 
-B, the behavior is different in the pnp and pp'p regimes, 
realized at positive and negative energies, respectively. In 
the pnp case, a crossover to FP oscillations is observed. 
In the pp'p case, discrete Landau levels give way to a con- 
tinuous spectrum in the region inside a parabola (dashed 
line) which marks the quasiclassical collapse threshold, 
Eq.®. 

We note that the behavior of dG/dVtg observed in 
the pnp region, in particular the half-period shift of FP 
fringes at relatively low B < IT, is not captured by local 
DOS, Eq.([Tni). As discussed in Ref.[ll|, this half-period 
shift results from FP interference due to Klein scatter- 
ing at pn interfaces. A proper model of this effect must 
account for ballistic conductance in the FP regime. 

The collapse observed in the density of states is related 
to deconfinement of classical orbits. The orbits can be 
analyzed as constant energy trajectories of the problem 



e = VF\/Pi +pI + U{x), Py =Py - eBx. (12) 

For parabolic U{x) = —ax^ the orbits with py = can 
be easily found in polar coordinates Px + ipy — |p|e*^: 
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FIG. 2: (a) Closed orbits for the Thomas-Fermi potential ob- 
tained from the density profile, Eq.((5]), with _B = 9, 7, 5, 3, 1 T 
and = 0. Long trajectories, extending far outside the gated 
region, do not contribute to SdH oscillations (see text). (b,c) 
Trajectories for the potential U(x) = —ax^ and Py = 0. Three 
types of trajectories are shown in momentum space (b) and 
position space (c): subcritical (red), critical (black), and su- 
percritical (blue). The saddle points in momentum space cor- 
respond to motion along straight lines x = it£, where the 
Lorentz force is balanced by the electric field. 



with po = vpe^B^ /2a (see Figl^jD). Only real, positive 
solutions should be retained; when e/sc > 1, the discrim- 
inant in Ea. (|13p is negative near sa 7r/2 and trajectories 
cannot close (blue curves in Figl2}D). 

The related orbits in position space can be found from 
the relation dy/dx = y/x = Py/Px, giving 



dx 



zLvp{py — eBx) 



v/(£ - U{x)f - vUpy - eBxf 



(14) 



For py = 0, integration is performed using the variable 

u = xye, 



y 



^ ± 



du 



V(M + £/ec-2)2 + 4(e/£,-l)' 



(15) 



where £ = VFeB/2a such that Ec = aP. The integrand 
changes its behavior at the critical energy Ec- For £ > £c, 
the integrand is real valued for all u and 



sinh 



(16) 



2 V^/£c - 1 

For £ < £c, real solutions are divided into two domains 
Q < u < 2 — e/£c — 2-y/l — £/£c (closed orbits) and u > 
2 — e/ec + 2\J\ — ejec (open orbits): 



cosh 



V{x) ~ yo 

I 



± 



2 - £/£, - x^ie 

2^/1 -e/ec 



(17) 



The red curves in Figl2]: correspond to the low energy 
regime, £ < £c, where orbits can either be closed (Landau 
levels) or open (trajectories for particle moving far from 
the barrier). At higher energies, £ > £c, all trajectories 
are open. The straight black lines correspond to the crit- 
ical orbits of Eq.([5]), where the Lorentz force and electric 



field are balanced. In addition to the two particular crit- 
ical trajectories shown, in the limit ejec — >■ 1 there is an 
entire family of critical trajectories which asymptotically 
approach these lines. 

Interestingly, unlike in the case of the potential ob- 
tained from the Thomas-Fermi model, where the clas- 
sical turning points move continuously to infinity as the 
transition is approached, trajectories in the parabolic po- 
tential are trapped between the critical separatrix lines. 
At very low energies, closed orbits are approximately cir- 
cular; as the energy increases towards £c, orbits become 
more and more elongated, until finally merging with the 
separatrix at £ = £c (see Figl^J . 

In summary, graphene devices with a barrier induced 
by a narrow top gate can be used to probe electronic 
states on the spatial scale of a few tens of nanometers. 
In our transport measurments, the SdH-type resonances 
arising from quantized states associated with closed or- 
bits are used to directly observe the competition between 
magnetic confinement and deconfinement due to elec- 
tric field. As a result of this competition, the discrete 
spectrum of Landau levels collapses when subjected to a 
strong external potential. Experimental observations are 
found to be in good agreement with theory. 
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